Mott— Hubbard transition vs. Anderson localization of correlated, disordered electrons 
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The phase diagram of correlated, disordered electrons is calculated within dynamical mean-field 
theory using the geometrically averaged ("typical") local density of states. Correlated metal, Mott 
insulator and Anderson insulator phases, as well as coexistence and crossover regimes are identified. 
The Mott and Anderson insulators are found to be continuously connected. 
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The properties of real materials are strongly influenced 
by the electronic interaction and randomness Q] . In par- 
ticular, Coulomb correlations and disorder are both driv- 
ing forces behind metal-insulator transitions (MITs) con- 
nected with the localization and delocalization of parti- 
cles. While the Mott-Hubbard MIT is caused by the 
electronic repulsion 0, the Anderson MIT is due to co- 
herent backscattering of non-interacting particles from 
randomly distributed impurities [jj. Furthermore, disor- 
der and interaction effects are known to compete in subtle 
ways 0, 0, • Several new aspects of this interplay will 
be discussed here. 

The theoretical investigation of disordered systems 
requires the use of probability distribution functions 
(PDFs) for the random quantities of interest. In phys- 
ical or statistical problems one is usually interested in 
"typical" values of random quantities which are mathe- 
matically given by the most probable value of the PDF 
0. In many cases the complete PDF is not known, i.e., 
only limited information about the system provided by 
certain averages (moments or cumulants) is available. In 
this situation it is of great importance to choose the most 
informative average of a random variable. For example, if 
the PDF of a random variable has a single peak and fast 
decaying tails this variable is usually well estimated by 
its first moment, known as the arithmetic average. How- 
ever, there are many examples, e.g., from astronomy, the 
physics of glasses or networks, economy, sociology, biol- 
ogy or geology, where the knowledge of the arithmetic 
average is insufficient since the PDF is so broad that its 
characterization requires infinitely many moments. Such 
systems are said to be non-self-averaging. One example 
is Anderson localization: when a disordered system is at 
the Anderson MIT [2j, most of the electronic quantities 
fluctuate strongly and the corresponding PDFs possess 
long tails Q ■ This is well illustrated by the local density 
of states (LDOS) of the system. The arithmetic mean 
of this random one-particle quantity does not resemble 
its typical value at all. In particular, it is non-critical 
at the Anderson transition y| and hence cannot help to 



detect the localization transition. By contrast, the geo- 
metric mean 0, H(ij , which gives a good approximation 
of the most probable ("typical ") value of the LDOS, 
vanishes at a critical strength of the disorder and hence 
rovides an explicit criterion for Anderson localization 
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A non-perturbative theoretical framework for the in- 
vestigation of correlated lattice electrons with a lo- 
cal interaction is given by dynamical mean-field the- 
ory (DMFT) HI 03. If in this approach the effect of 
local disorder is taken into account through the arith- 
metic mean of the LDOS one obtains, in the absence 
of interactions, the well known coherent potential ap- 
proximation |l7j |. which does not describe the physics of 
Anderson localization. To overcome this deficiency Do- 
brosavljevic and Kotliar formulated a variant of the 
DMFT where the geometrically averaged LDOS is com- 
puted from the solutions of the self-consistent stochastic 
DMFT equations. Employing a slave-boson mean-field 
theory as impurity solver they investigated the disorder- 
driven MIT for infinitely strong repulsion off half-filling. 
Subsequently, Dobrosavljevic et al. [12( incorporated the 
geometrically averaged LDOS into the self-consistency 
cycle and thereby derived a mean-field theory of Ander- 
son localization which reproduces many of the expected 
features of the disorder-driven MIT for non-interacting 
electrons. This scheme uses only one-particle quantities 
and is therefore easily incorporated into the DMFT for 
disordered electrons in the presence of phonons 0], or 
Coulomb correlations. 

In this Letter we employ the DMFT with the typical 
LDOS to determine the non-magnetic ground state phase 
diagram of the disordered Hubbard model at half-filling 
for arbitrary interaction and disorder strengths. Thereby 
the Mott-Hubbard and Anderson MITs are investigated 
on equal footing. The system is described by a single- 
orbital Anderson-Hubbard model 
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where t > is the amplitude for hopping between near- 
est neighbors, U is the on-site repulsion, ni a — a\ a a ia is 
the local electron number operator, a,i a is the an- 

nihilation (creation) operator of an electron with spin c, 
and the local ionic energies e, are independent random 
variables. In the following we assume a continuous prob- 
ability distribution for a, i.e., Vfe) = 6(A/2 - |ei|)/A, 
with G as the step function. The parameter A is a mea- 
sure of the disorder strength. 

This model is solved within DMFT by mapping it 0] 
onto an ensemble of effective single-impurity Anderson 
Hamiltonians with different e^: 

-ffsiAM = >Xei - /^atyffljg + Uriifnn (2) 

ko" k<x 

Here /i = —U/2 is the chemical potential correspond- 
ing to a half-filled band, and Vk and ek are the hy- 
bridization matrix element and the dispersion relation of 
the auxiliary bath fermions Ck CT , respectively. For each 
ionic energy we calculate the local Green function 
G(u>, ej), from which we obtain the geometrically aver- 
aged LDOS p gcom (w) = exp [(In Pi (co))] HH01G3, where 
Pi(uj) = -lmG(uj,e l )/TT, and (Oj) = / deiP{e i )0(e i ) is 
the arithmetic mean of Oi |2fl( ■ The lattice Green func- 
tion is given by the corresponding Hilbert transform as 
G(oj) = J du>' /Cgeom(w) / (to — lo'). The local self-energy 
is determined from the k-integrated Dyson equa- 
tion £(w) = lo — t](uj) — G _1 (w) where the hybridization 
function t](lo) is defined as 77(0;) = J2k l^k| 2 / (w — £k)- 
The self-consistent DMFT equations are closed through 
the Hilbert transform G(lo) = J deN (e)/ [u-e - T,(co)], 
which relates the local Green function for a given lat- 
tice to the self-energy; here Nq(c) is the non-interacting 
DOS. We note that this approach describes only the ex- 
tended states since the localized part of the spectrum, 
given by the isolated poles of G(lo), is not included. For 
this reason p ge0 m{oj) is not normalized to unity. 

The Anderson-Hubbard model is solved for a semi- 
elliptic DOS, N (e) = 2y/D 2 ~ e 2 /irD 2 , with bandwidth 
W = 2D; in the following we set W = 1. For this 
DOS a simple algebraic relation between the local Green 
function G{lo) and the hybridization function T](u>) = 
D 2 G{lo)/A holds [3. The DMFT equations are solved at 
zero temperature by the numerical renormalization group 
(NRG) technique |2l|,|22j which allows us to calculate the 
geometric average of the LDOS in each iteration loop. 

The main result of this Letter is the ground state phase 
diagram of the Anderson-Hubbard model at half-filling 
shown in Fig. ^ Two different phase transitions are 
found to take place: a Mott-Hubbard MIT for weak dis- 
order A, and an Anderson MIT for weak interaction U. 
The two insulating phases surround the correlated, dis- 
ordered metallic phase. The properties of these phases, 
and the transitions between them, will now be discussed. 
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FIG. 1: Non-magnetic ground state phase diagram of the 
Anderson-Hubbard model at half-filling as calculated by 
DMFT with the typical local density of states. 
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FIG. 2: Local density of states (LDOS) as a function of disor- 
der A for various values of the interaction U. Solid (dashed) 
curves correspond to the geometrically (arithmetically) aver- 
aged LDOS. 



(i) Metallic phase: The correlated disordered metallic 
phase is characterized by a non-zero value of /o ge om(0), 
the spectral density at the Fermi level (lo = 0). With- 
out disorder DMFT predicts this quantity is to be given 
by the bare DOS iVo(O), as expressed by the Luttinger 
theorem |2^ | . This means that Landau quasiparticles are 
well-defined at the Fermi level. The situation changes 
dramatically when randomness is introduced, since a sub- 
tle competition between disorder and electron interaction 
arises. Increasing disorder at fixed U reduces p ge om(0) 
and thereby decreases the metallicity as shown in the 
upper panel of Fig. [5J The opposite behavior is found 
when the interaction is increased at fixed A (see Fig. |3 
for A = 1), i.e., the metallicity improves in this case. 
In the strongly interacting metallic regime the value of 
Pgeom(O) is restored, reaching again its maximal value 
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FIG. 3: Geometrically averaged LDOS as a function of inter- 
action U for different disorder strengths A. Solid (dashed) 
curves with closed (open) symbols are obtained with an ini- 
tial metallic (insulating) hybridization function. Triangles: 
A = 1; dots: A = 2.5. Left inset: LDOS with Mott gap at 
U — 3 for different disorder strengths A. Right inset: Inte- 
grated LDOS iVext as a function of A at U = 3 and 5. 
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FIG. 4: LDOS for U = 1.25 (left column) and U = 1.75 (right 
column) for different disorder strengths A. Solid (dashed) 
curves correspond to the geometrically (arithmetically) aver- 
aged LDOS. 

Nq(0). Physically this means that in the metallic phase 
sufficiently strong interactions protect the quasiparticles 
from their decay due to impurity scattering. For weak 
disorder this effect of the interaction is essentially inde- 
pendent of the choice of the LDOS. 

(ii) Mott-Hubbard MIT: For weak to intermediate dis- 
order there is a sharp transition at a critical value of U 
between a correlated metal and a gapped Mott insula- 
tor. We find two transition lines depending on whether 
the MIT is approached from the metallic side [A^ H (U), 
full dots in Fig. P or from the insulating side [A^[ H (U), 
open dots in Fig. [U . This is very similar to the case with- 
out disorder 0, 122. Eij ; the hysteresis is clearly seen in 
Fig. for A = 1. The Af x H (U) and A^ H (U) curves in 



Fig. ^have positive slope. This is a consequence of the 
disorder-induced increase of spectral weight at the Fermi 
level (see Fig. QJ which in turn requires a stronger inter- 
action to open the correlation gap. In the Mott insulating 
phase close to the hysteretic region an increase of disor- 
der will therefore drive the system back into the metallic 
phase. The corresponding abrupt rise of p gcom (0) is seen 
in the lower panel of Fig. [3 and the right column of 
Fig. 0J In this case the disorder protects the metal from 
becoming a Mott insulator. 

Around A w 1.8 the A* m (C7) and A^ H (U) curves ter- 
minate at a single critical point, cf. Fig. ^ At stronger 
disorder (A > 1.8) only a smooth crossover from a metal 
to an insulator takes place. This is clearly illustrated 
by the U dependence of /o ge om(0) shown in Fig. [31 for 
A = 2.5. In this parameter regime the Luttinger theo- 
rem is not obeyed for any U. In the crossover regime, 
marked by the hatched area in Fig. ^ Pgeom(O) vanishes 
gradually, so that the metallic and insulating phases can 
no longer be distinguished rigorously |25j . 

Qualitatively, we find again that the Mott-Hubbard 
MIT and the crossover region do not depend much on 
the choice of the average of the LDOS. We also note the 
similarity between the Mott-Hubbard MIT scenario dis- 
cussed here and the one for the system without disorder 
at finite temperatures (UEHIill, especially the presence 
of a coexistence region with hysteresis. However, while in 
the non-disordered case the interaction needed to trigger 
the Mott-Hubbard MIT decreases with increasing tem- 
perature, the opposite holds in the disordered case. 

(Hi) Anderson MIT: In Fig. ^the metallic phase and 
the crossover regime are seen to lie next to an Anderson 
insulator phase where the LDOS of the extended states 
vanishes completely. The critical disorder Af(U) cor- 
responding to the Anderson MIT is a non-monotonous 
function of the interaction: it increases in the metallic 
regime and decreases in the crossover regime. Where 
A^(U) has a positive slope an increase of the interac- 
tion turns the Anderson insulator into a correlated metal. 
This is illustrated in Fig. 02 for A = 2.5: at U/W » 0.7 
a transition from a localized to a metallic phase occurs, 
i.e., the spectral weight at the Fermi level becomes fi- 
nite. In this case the electronic correlations impede the 
localization of quasiparticles due to impurity scattering. 

Fig. [21 shows that the Anderson MIT is continuous. 
In the critical regime /o geom (0) ~ [A^ (U) — A]' 3 for 
U = const. In the crossover regime we find a critical 
exponent (3—1 (see the case U = 1.75 in lower panel 
of Fig. EJ; elsewhere /3 ^ 1. However, since it is diffi- 
cult to determine (3 with high accuracy we cannot rule 
out a very narrow critical regime with (3 = 1. It should 
be stressed that an Anderson transition with vanishing 
Pgeom(O) at finite A = A^(U) can only be detected in 
DMFT when the geometrically averaged LDOS is used 
(solid lines in Fig. With arithmetic averaging one 
finds a nonvanishing LDOS at any finite A (dashed lines 
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in Fig. 01. 

(iv) Mott and Anderson insulators: The Mott insula- 
tor with a correlation gap is rigorously defined only for 
A = 0, and the gapless Anderson insulator only for U = 
and A > A^(0). In the presence of interaction and dis- 
order this distinction can no longer be made. However, 
as long as the LDOS shows the characteristic Hubbard 
subbands (see left inset in Fig. one may refer to a 
disordered Mott insulator. With increasing A the spec- 
tral weight of the Hubbard subbands vanishes (see right 
inset in Fig. [3J and the system becomes a correlated An- 
derson insulator. The border between these two types 
of insulators is marked by a dashed line in Fig. ^ The 
results obtained here within DMFT prove that the para- 
magnetic Mott and Anderson insulators are continuously 
connected. Hence, by changing U and A it is possible to 
move from one type of insulator to the other without 
crossing the metallic phase. 

In conclusion, using DMFT with the geometrically av- 
eraged (typical) LDOS we computed the non-magnetic 
ground state phase diagram of the Anderson-Hubbard 
model at half-filling for arbitrary interaction and disor- 
der strengths. In particular, we determined the position 
of the Mott-Hubbard metal-insulator and Anderson lo- 
calization transitions. The presence of disorder increases 
the critical interaction for the Mott-Hubbard MIT, and 
turns the sharp transition (with hysteresis) into a smooth 
but rapid crossover. On the other hand, the critical disor- 
der strength for Anderson localization increases for weak 
interaction and is suppressed by strong interactions. The 
paramagnetic Mott and Anderson insulators are contin- 
uously connected. The specific predictions of our the- 
ory not only apply to disordered solids but also to cold 
fermionic atoms in optical lattices (2?| ■ In the latter case 
a precise control of system parameters appears to be pos- 
sible which, in principle, allows one to explore all parts 
of the phase diagram. 
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